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ABSOLUTELY MINIMISING GENERALISED SOLUTIONS TO 
THE EQUATIONS OF VECTORIAL CALCULUS OF 
VARIATIONS IN 

NIKOS KATZOURAKIS 


Abstract. Consider the supremal functional 

(1) _Eoo(w,A) := ACn, 

applied to maps u : Q C R —> R^, A > 1. Under certain assumptions 

on Jf, we prove for any given boundary data the existence of a map which is: 

i) a vectorial Absolute Minimiser of (1) in the sense of Aronsson, 

ii) a generalised solution to the ODE system associated to (1) as the analogue 
of the Euler-Lagrange equations, 

iii) a limit of minimisers of the respective functionals as p —>■ oo for any 
g > 1 in the strong W^’'^ topology &; 

iv) partially on H off an exceptional compact nowhere dense set. 

Our method is based on approximations and stable a priori partial regular¬ 
ity estimates. For item ii) we utilise the recently proposed by the author notion 
of D-solutions in order to characterise the limit as a generalised solution. Our 
results are motivated from and apply to Data Assimilation in Meteorology. 
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1. Introduction 

Calculus of Variations in has a long history and was pioneered by Aronsson 
in the 1960s [Al]-[A5]. In the vector case and in one spatial dimension, the basic 
object of study is the functional 

(1.1) E^{u,A) ■= u-.ncR — Acn. 
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Here u G iV > 1, H is an open interval, A is measurable and AC G 

(n X X K^) is a function which we call Lagrangian and whose arguments will 
be denoted by (x, rj, P). Aronsson who studied the case N = 1 was the first to note 
the locality problems associated to this functional. By introducing the appropriate 
minimality notion in among other things proved the equivalence between the 
so-called Absolute Minimisers and solutions of the analogue of the Euler-Lagrange 
equation which is associated to the functional under (7^ smoothness hypotheses. 
The minimality notion of Aronsson adapted to the vector case of (1.1) is 

(1.2) E^{u,n') < E^{u + ct>,n'), VH' (£ ((>G W'o’~(H',K^). 

The higher dimensional scalar analogue when u : C K" —>■ K is a real function 
has also attracted considerable attention by the community and by now there is 
a vast literature, for instance see Crandall [C], Barron-Evans-Jensen [BEJ], and 
for a pedagogical introduction see [KO] and references therein. In particular, the 
Crandall-Ishii-Lions theory of Viscosity solutions proved to be an indispensable tool 
in order to study the equations in L°° which are non-divergence, highly nonlinear 
and degenerate. Even in the simplest case where the Lagrangian is the Euclidean 
norm, i.e. A^'{P) = \P\'^, in general the solutions are non-smooth and the corre¬ 
sponding PDE which is called oo-Laplacian for smooth functions reads 

n 

(1.3) AooM := Du(8)Du:D^u = ^ DiuDjuD^jU = 0. 

jj'=i 

However, until the early 2010s, the theory was essentially restricted to the scalar 
case N = 1. A most notable exception is the early vectorial contributions of Barron- 
Jensen-Wang [BJWl, BJW2]. Therein the authors among other far-reaching results 
studied the weak* lower semicontinuity of general supremal functionals and proved 
under certain assumptions the existence of Absolute Minimisers in the “rank-1” 
cases, i.e. when either n = 1 or V = 1. 

In a series of recent papers [KI]-[K9], the author has initiated the systematic 
study of the vector-valued case, which except for its intrinsic mathematical inter¬ 
est, appears to be important for many real-world applications (see also the joint 
contributions with Abugirda, Croce, Pisante and Pryer [AK, CKP, KP, KP2]). In 
particular, the complete PDE system arising in L°° was derived and studied in 
[Kl]. The results in the aforementioned papers include in particular the study of 
the analytic properties of classical solutions to the fundamental equations and their 
connection to the supremal functional. In the case of 

(1.4) Eoo(u,A) = — >R^,Acn, 

(where |Du| denotes the Euclidean norm of the gradient on K-^^”), the respective 
oo-Laplace system written for smooth maps is 

(1.5) AooM := (g) Du -|- iDuppu]-*-®: D^u = 0. 

In (1.5), [Du(a;)]"*~ denotes the orthogonal projection on the orthogonal complement 
of the range of the linear map Du(a;) : K" —and in index form reads 

N n 

{piUaPijUp + |Dup[Du]^^(5y^ DyM/3 = 0, a = 1,..., A, 

/3=1 i ,i=l 


[Du]-L := Proj^(D„)_L. 
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An extra difficulty of (1.5) which is not present in the scalar case of (1.3) is that the 
coefficients may be discontinuous along interfaces even for C°° solutions because the 
term involving [Du]-*- measures the dimension of the tangent space of u{n) C 
(see [Kl, K2] and the numerical experiments in [KP]). This is a general vectorial 
phenomenon studied in some detail in [K.3]. The appropriate minimality notion 
allowing to connect (1.5) to the functional (1.4) has been established in [K4]. It 
is a remarkable fact that when the rank of the gradient is greater than one, i.e. 
for maps u : C K" —y such that min{n, A} > 2, Absolute Minimimality 

for (1.4) is neither necessary nor sufficient for solvability of (1.5) and the correct 
notion of oo-Minimal maps is intrinsically different (see [K4]). In the case of 
the supremal functional (1.1), the associated equations written for smooth maps 
u : n C M — y read 


( 1 . 6 ) 

where 


(1.7) 


(•, M, Du,D^u) = 0, on n, 


Tao{x,ri,P,X) 


Af’p{x, rj, P) ® Afp{x, rj, P) 

+ Af'{x,ri,P)[Af’p{x,ri,P)]-^Af'pp{x,r],P) 


X 


+ (^A!fn{x,r],P) ■ P + Afx{x,ri,P)^^p{x,ri,P) 
+ -^{x, V, P) [^pix, r], P)] (Afp^(x, 77, P)P 
+ .^Px{x,r],P) - Af’rf{x,r],P)y 


In (1.7), the notation of subscripts denotes derivatives with respect to the respective 
variables and [Jifp{x,ri, P)] is the orthogonal projection 

(1.8) [Af’p{x,r],P)]'^ := 7- sgn(.ifp(x,r 7 ,P)) (g)sgn(ifp(a;,77,P)). 

The system (1.6)-(1.8) is a 2nd order ODE system which is quasilinear, non¬ 
divergence, non-monotone and with discontinuous coefficients. Even in the scalar 
case of A = 1 in which Poo simplifies to 

Pao{x,r],P,X) = {^p{x,r],P)fx + (j!f^{x,r], P)P + J!f,,{x,r], P)^Jf'p{x,r], P) 

it is known since the work of Aronsson that in general does not have solutions 
any more regular than at best C'^(D,IR'^) and their “weak” interpretation is an 
issue. Let us also note that, also inspired by Aronsson’s work, Sheffield-Smart [SS] 
made a vectorial breakthrough relevant to (1.5) and (1.4) which was simultaneous to 
[Kl]. They studied smooth vector-valued optimal Lipschitz extensions of functions, 
deriving a different more singular version of oo-Laplacian than (1.5), corresponding 
to (1.4) but when the matrix norm of Du is the nonsmooth operator norm on 

In this paper we study the functional (1.1), the associated nonlinear system (1.6)- 
(1.8) and their connection. Our main result establishes for any given endpoint data 
on D the existence of a vectorial Absolute Minimiser u°° € IEj’°°(D,K^) of (1.1) 
which also is a generalised solution to (1.6)-(1.8) in a certain new sense to be made 
precise below. We moreover glean extra information about u°°; it a partially 
map off a singular compact set and also is a limit of minimisers of the respective 
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L™ functionals 

(1.9) Emiu^n) := u-.ncR—> , 

Jn 

in the strong W^’'^ topology as m —)■ oo for any q > 1. Our results have been 
motivated from and apply to variational Data Assimilation in the form used in 
the Earth Sciences and in particular Weather Forecasting. Below we discuss the 
essential idea of our new notion of generalised solution, our assumptions, state our 
main result and also draw connections to Data Assimilation. 

Motivated in part by the nonlinear systems arising in in the very recent 
paper [K8] the author proposed a new theory of generalised solutions which applies 
to fully nonlinear PDE systems. In addition, this theory allows to interpret merely 
measurable general mappings u : C K" — y as solutions of PDE systems 

which may even be defined by discontinuous nonlinearities and can be of any order. 
Our approach is duality-free and bypasses the insufficiency of the standard duality 
ideas to apply to even linear non-divergence equations with rough coefficients. The 
standing idea of the use of integration-by-parts in order to pass derivatives to test 
functions is replaced by a probabilistic description of the limiting behaviour of the 
difference quotients. This builds on the use of Young measures valued into compact 
tori, which is the compactification of the space wherein the derivatives are valued. 
Background material on Young measures can be found e.g. in [EG, CFV, V, EL, 
P, E, M], but for the convenience of the reader we recall herein the rudimentary 
properties we actually utilise. 

The essential idea of our new notion of solution for the case needed in this paper 
can be briefly motivated as follows. Assume that u : td CR — R^ is a strong a.e. 
solution of the system 

(1.10) E[-,u,Du,D'^u) = 0, on D, 

in Wj^’“(D, M^). We need a notion of solution which makes sense even if u is 
merely W[Q^(r2,K'^). To this end we rewrite the above statement that u is a 
strong solution in the following unconventional fashion 


( 1 . 11 ) 


sup 

AGsupp((5j;,2„(„, 


E(x, u{x), Dm(x), X 


= 0 , 


a.e. X G fl. 


That is, we change from the classical viewpoint that the 2nd derivative is a map 
D^u : n C M — R^ valued in R^ to that it is a probability valued map given by 
the Dirac mass at TEw. 


5d 2„ : n C K - > ^{R^), Xi -(5D2u(a;)- 

Obviously, “supp” denotes the support of the probability measure. Similarly, if D^’^ 
stands for the difference quotient operator, it can be shown that we may rewrite 
the definition of D^u as 


(1.12) (5i3i,hp)„ —^ (5 d2u, as/i —> 0. 

The weak* convergence above is meant in the so-called Young measures into R^, 
that is the probability-valued maps 4 ? : D C M —>• !^{R^) which are weakly* 
measurable (see the next section for details and Lemma 2.3). The rationale behind 
the reformulation (1.11)-(1.12) is that we may allow more general probability-valued 
maps arising as “diffuse” 2nd derivatives for maps classically differentiable only 
once. The latter of course may no longer be the concentration measure-valued 
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maps This is indeed possible if we augment and replace it by its 1- 

point compactification U {oo}. By considering Young measures valued into 
spheres, we obtain the necessary compactness and the maps (( 5 Di,feDii)ft. 5 ^o always 
have subsequential weak* limits 

(5^)1, u, as hi —^ 0 

in the space of sphere-valued Young measures fl C M —U {oo}) (even 
when u is merely once differentiable). The above ideas motivate the notion of V- 
solutions and diffuse derivatives in the special case of W^’°° solutions to (1-10) 
(see Definitions 3.1 and 3.2) and will be taken as principal in this work. 

As a first application of this new approach, in the paper [K8] among other 
things we proved existence of I?-solutions to the Dirichlet problem for (1.5) when 
n = N. Further results in the context of X>-solutions have been established in 
[K9, KIO, Kll], including certain uniqueness assertions. Herein we focus on (1.1) 
and (1.6)-(1.8). This is a non-trivial task even in the ID case. In fact, it is not 
possible to work in the generality of (1.1), (1.6)-(1.8) without structural conditions 
on The most important restriction is that the Lagrangian has to be radial in 
P. This means that .if can be written as 

(1.13) ^{x,T],P) = (^x,r],^\P - 'r{x,r])\‘^^ . 


for some mappings Jif : O x x [0,oo) —> M and Y : H x R-^ —> R^. This 
condition is justified by the results of [K3] since as we proved therein it is both 
necessary and sufficient for the ODE system to be degenerate elliptic. The extra 
bonus is that then the coefficients of (1.7) match and become continuous. This 
is a special occurrence due to the ID nature of the problem and can not happen 
when n > 2. Anyhow, under the assumption (1.13), (1.7) after a certain rescaling 
becomes 


(1.14) 


Poo{x,T],P,X) = \P - y{x,r])\ JTp [x,r],-\P - fix,?])] I X - %,{x,-ri)P 


-%,{x,r]) + {P -f'{x,r]))jifp[x,r],-\P -•r{x,r]) 


JYf {x,rj,f-\P -f'{x,r])r ) -|- P ■ JYp [ x,r],l-\P - y{x,r]) 


- jif;{x,v,-\p-y{x,p)n \p-y{x,r,)n 




- [P -P{x, rj)) ® [P -f {x, rj)) 

- {P- y{x,ri)Y'Pp{x,-q)jep (^,'q,]^\P - y{x,Ti)'f^ 


The naturalness of our structural assumption (1.13) is also justified by Lagrangian 
models arising in variational Data assimilation which we describe briefly after the 
statement of our main result. 
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Theorem 1.1. Let C K 6e a bounded open interval and let also 
-W : nx]R^x[0,oo)—^[l,oo), ^ : H x 


be given maps with N We suppose that 


(1.15) < 

and also 

(1.16) 


is up to the boundary, 

C{\r]\) > J^p{x,r],p) > Co, 
2Jifpp{x,r],p)p + Jfp{x,r],p) > Cq, 

\jf,,{x,7],p)\ + \ji^p{x,r],p)\ < Cd??!))! + p), 
J^pp{x,r],p)\ + \jifpp{x,r],p)\ + \j^p^{x,r],p)\ < C{\r]\){l + p^), 


'y is up to the boundary, 

< {l/co)(l + IvH, 


for some constants co,Q! G (0,1), some M gN, some positive continuous increasing 
function C G C^([0,oo)) and all {x,r],p) G LI x x [0, oo). Then, for any affine 
map b : K — > R^, there exists a map u°° G K^^) with the following 

properties: 

(1) is a vectorial Absolute Minimiser of the functional 

(1.17) Eao{u,A) = esssvLpJ^[x,u{x),\\T)u{x) — 'y{x,u{x))'\^], 

xgA V 2 ' J 


that is it satisfies (1.2). 

(2) u°° is a D-solution (see Definitions 3.1 and 3.2) of the system 

(1.18) J^oo(-,m,Dm,D^u) = 0, on LI, 
where Too is given by (1.14). 

(3) u°° is a subsequential limit as m —t' oo in the strong K^) topology 

of minimisers {u ^}^^2 of the functionals 

(1.19) E,^{u,A) = f jy 

JA 

where each minimises over the respective space K^), for any 

q G [l,cx)). 

(4) There is an open subset Ltoo L= LI such that u°° G C^(Llao,R^). Moreover, 

Ll\Ll^= 5({Du°° = r(-,u“)}) 

and hence LI \ Ll^o is compact and nowhere dense in LI. 


1 I 12 

x,u{x), - DM(a;) — y{x,u{x))\ 


dx. 


Item (4) above is a partial regularity assertion which differs from more classical 
results in that the singular set LI \ Liao is a relatively compact nowhere dense set (a 
topological boundary) but not necessarily a Lebesgue nullset. This is a new type of 
partial regularity which seems to arise in L°°. Item (3) indicates the fashion in which 
u°° is obtained, namely via the well-established method of L™ approximations of 
problems as m ^ oo, but also includes a non-trivial fact, the strong convergence 
of the L"* minimisers u™ together with their first derivatives to u°°. Note that 
Theorem 1.1 above does not state that (1) implies (2), but instead that there is an 
object u°° which satisfies both. In order to obtain solely (1), the hypotheses (1.15)- 
(1.16) can be relaxed substantially (accordingly, see the paper [AK]), but since 
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herein we are interested in the satisfaction of the equations as well we do not tackle 
this problem separately. Finally, due to the dependence of ^ on the lower order 
terms (a;, u(x)), the Absolutely Minimising I?-solution u°° is not in general unique, 
as shown by the example ^{x,P) = sin^ x |Pp of Yu [\ ] even when n = N = 1 
and P = 0. Uniqueness is a most delicate question already in the scalar case (see 
[.JWY, MWZ]). Let us also recall that (1) above has been obtained in [BJWl] but 
under different hypotheses on ^ which in particular require J!f{x,r],0) = 0 and 
"U = 0, a fact incompatible with “additive” Lagrangians like those arising in Data 
Assimilation which we describe right next. 

The motivation to study the present ID vectorial L°° variational problem comes 
from Data Assimilation models arising in the Earth Sciences and especially in Me¬ 
teorology. More precisely, following the terminology of [B], we are inspired by 
a continuous time generalisation of what is known as weakly constrained four¬ 
dimensional variational assimilation (4D-Var) in the geosciences. For more details 
we refer e.g. to [AJSV, BP, CT, De, FS, LDT, PVT, RCPTV, RJ, TC, Y]. 

Let us describe briefly the model in pure mathematical terms. Let Y : U x 
be a time-dependent vector field describing the law of motion of a 
body moving along a trajectory defined by the solution u : U C K —of Du = 
'V(-,u) (e.g. Newtonian forces, finite-dimensional Galerkin approximation of the 
Euler equations, etc.). Let also k : D CR — R^ be some partial “measurements” 
in continuous time along the trajectory and K : R^ — > be a submersion 

which corresponds to some component of the trajectory we are able to measure, for 
example some projection. Then, we wish to find a u which should satisfy the law 
of motion and also be compatible with the measurements along the trajectory: 

Du{t) = y(t,u{t)) & K{u{t)) = k{t), t G n. 

Lfowever, this problem is in general overdetermined (due to errors in the measure¬ 
ments, etc.) since we impose a pointwise constraint to the solution of the system. 
In standard variational Data Assimilation (see [B, BS]), one instead seeks for ap¬ 
proximate solutions by minimising the “error” integral functional Ei given by (1.9) 
for TO = 1 and with Lagrangian given by 

(1.20) ^{x,p,P) := 1 + ^\k{x) - K{p)\^ + ^|P-r(T,77)|' 

which describes the “error”. But if instead we choose to use the respective supremal 
functional (1.1) with .if as in (1.20), large “spikes” of the deviation from the actual 
solution with small area are from the outset excluded. For the Lagrangian (1.20), 
the equations (1.6)-(1.14) arising in Data Assimilation read 


- [Du-P(-,u)]^ 

(1.21) . ((iF(u) - k^K^iu) - (iru)^r^(-,u)) - [x^iu) : {K{u) - k) 

® Du -I- {K{u) — kfj ■ kx (Du — P (•, u)) = 0 


Our main result applies in particular to (1.20)-(1.21). Although the equations 
are more complicated than the respective Euler-Lagrange equations, evidence 
obtained in [BK] suggests that they provide more accurate models. 
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2. Basics and Young measures into spheres 

Our basic notation is either standard (as e.g. in [E2, EG]) or self-explanatory. 
Eor example, the Lebesgue measure on K will be denoted by | • |, the characteristic 
function of the set Ahy XAi the standard Sobolev and spaces of maps u : O C 
K —^ by LP(n,]R^), etc. We will also follow the standard 

practice that while deriving estimates, universal constants may change from line to 
line but will be denoted by the same letter. Y S N will always be the dimension of 
the range of our candidate solutions u : O C M —Unless indicated otherwise, 
Greek indices a,/3,7,... will run in {1,..., Y} and the summation convention will 
be employed in products of repeated indices. The standard basis on will be 
denoted by {e^,..., e^} and hence for the map u with components Ua we will write 
u{x) = Ua(x)e°‘. The norm symbol j • j will always indicate the Euclidean one and 
the respective inner product will be denoted by Given ^ S K^, we define for 
later use the following orthogonal projections of 

(2-1) := sgn(0 Gsgn(f), [^j-^ := I - sgn(^) (g) sgn(^). 

Here “sgn” stands for the sign function: sgn(^) := ^/j^j when ^ ^ 0 and sgn(O) := 0. 

Let now Y C M be a (Lebesgue) measurable set and consider the Alexandroff 
1-point compactification of the space 

:= K^U{oo}. 

Its topology will be the standard one which makes it isometric to the Y-sphere 
(via the stereographic projection which identifies {oo} with the north pole). The 
space will be considered equipped with the metric topology induced by the 
embedding into its compactification but balls, distances, etc. will be taken with 
respect to its usual metric structure. 

Definition 2.1 (Young Measures into the 1-point compactification of K^). The 
space of Young Measures E C R — i- R^ is denoted by (Y, R^) and is the set of 
probability-valued maps 

raebx^ -dix) e 

which are measurable in the following sense: for any fixed open set U C R^, the 
function Y 9 x i —> [i?(x)](Y) € K is measurable. This is called weak* measurability. 

The set (Y, R^) is a subset of the unit sphere of the space Y“, (Y, Ai (R^)). 
This Banach space consists of weakly* measurable maps valued in the signed Radon 
measures: Y 9 x i—i?(x) G The norm of the space is 

•= ess sup y(x)||(]R^) 

^ xGE 

where “|| • ||(K^)” is the total variation. For more details about this and relevant 
spaces we refer e.g. to [EL] (and references therein). Hence, the Young Measures 
are the subset of the unit sphere which consists of probability-valued weakly* mea¬ 
surable maps. It can be shown (see e.g. [EL]) that Yjjjl (Y, (K'^)) is the dual 

space of the space of measurable maps valued in the (separable) space C°(R^) 
of real continuous functions over R^, in the standard Bochner sense: 

{L\E,C°{R^))y = Y~ (Y,M(1^)). 
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The elements of this space are certain Caratheodory functions $ : if x ^ K 

endowed with the norm 

The space L^(if,C'°(K'^)) is separable and the duality pairing 

(•,•) : X L\E,C°{W^)) 

is given by 

:= [ [ ^{x,X)d[i){x)]{X)dx. 

JE Jm"' 

The unit ball of (if, Af is sequentially weakly* compact. Hence, for 

any bounded sequence C L“. (if, Ai , there is a limit map d and a 

subsequence of m’s along which as m —>■ oo. 

Remark 2.2 (Properties of ^^(if,K'^)). The set of Young measures is convex and 
by the compactness of R^, it can be proved that it is sequentially weakly* compact 
in (if, Ai(M^)) (see e.g. [FG, CFV]). This property is essential in our setting. 
Moreover, the space of measurable maps u : if C K — y R^ can be (nonlinearly) 
embedded into i^(if,K^) by v i—(5„, {Sy){x) := Sy(^x}, x G E. 

The following lemma is a small variant of a standard result about Young measures 
but it plays an important role in our setting (for the proof see [FG, GFV, V, K8]). 

Lemma 2.3. Let v'^,v°° : if C K — R^ be measurable maps, m G N. Then, up 
to the passage to subsequences, the following equivalence holds true as m ^ oo: 

ym — yv°°,a.e.onE ^ Syr,. ^ 6y^, in {E,R^). 

3. P-SOLUTIONS FOR FULLY NONLINEAR SYSTEMS 

Now we give the main definitions of our notion of solution only in the special 
case which is needed in this paper. For the general case and applications we refer 
to [K8]-[K11] and [GKP, KP, KP2]. 

Definition 3.1 (Diffuse 2nd derivatives). Suppose that u : LI C R — R^ is in 
Wjq“(D,K^). For any /i ^ 0, we consider the difference quotients of the derivative 

D1’'‘Du = ^ (du(- + h)- Du) : D C K — y R^ 

and Du is understood as being extended by zero on M \ 11. We define the diffuse 
2nd derivatives of u as the subsequential limits 'D'^u of ^di>d« in the space of 
Young measures from D into R^ along infinitesimal sequences Q R \ {0}: 

in ^^(H, R^), as i—y oo. 

The weak* compactness of ^^(ll,R'^) implies that every u G kF[g“(D,R^) pos¬ 
sesses diffuse 2nd derivatives, in particular at least one for every choice of (hi)^. 
For our notion of generalised solution, let us first introduce the following notation: 
if 1 ? is a probability measure on R'^, we define its reduced support as 

supp*(i?) := supp(i?) \ {oo} C R^ . 
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Definition 3.2 (I?-solutions of 2nd order ODE systems). Let 

(3.1) T ■. (O C K) X X X 


be a Borel measurable map with O an open set. Consider the ODE system 


(3.2) J^(-, u, Dm, D^m) = 0, on O. 

We say that a map m : O C M — >• in R-^) is a I?-solution of (3.2) 

when for any diffuse 2nd derivative 'D'^u € ‘S/ (O, we have 


(3.3) 


sup 

X Gsupp^ {T?^u{x)) 


T 


(^x, u{x),Du{x) 


X 


= 0 , 


a.e. X G fl. 


In general diffuse derivatives may not be unique for non-differentiable maps. 
Moreover, (3.3) is trivially satisfied at certain points at which it may happen that 
'D‘^u{x) = (5{oo} and hence supp*(I?^M(a:)) = 0 (see also the examples in [K8, Kll, 
KIO]). It is an immediate consequence of Lemma 2.3 that diffuse derivatives are 
compatible with classical derivatives, in the sense that if u is twice differentiable a.e. 
on fl, then the diffuse 2nd derivative 'D'^u is unique in the sense that V^u = 5x^2^ 
a.e. on fl. The converse is true as well if D^m is interpreted in the sense of Ambrosio- 
Maly as “derivative in measure” (see [AM, K8]). As a direct consequence we have 
that I?-solutions are compatible with a.e. twice differentiable strong solutions. 


4. The fundamental equations arising in L°° 

In this section we formally derive the fundamental equations (1.6)-(1.8) and 
in particular (1.14) associated to L°° variational problems for (1.1). The formal 
derivation of (1.6)-(1.8) has been performed in [Kl], but we include it here because 
it provides insights of the method of proof which makes the foregoing calculations 
rigorous. We obtain the L°° equations in the limit of the Euler-Lagrange equations 
related to the L"* integral functional (1.9) as m —oo. Here we suppose that m >2. 
The Euler-Lagrange equation of (1.9) is the ODE system 

(4.1) D(.if'"-i(-,M,DM).J^p(-,M,DM)) = .if'"-i(-,M,DM).if^(-,M,DM). 

By distributing derivatives and normalising, (4.1) gives 

(4.2) 

D(.if (•, M, DM)).ifp(-, M, Dm) -|- ———^ ^D(.ifp(-, m. Dm)) — m, Dm)^ = 0. 

Then, by employing (2.1) applied to ^ = .j2p(-,m,Dm) we expand the term in the 
bracket of (4.2) and obtain 

D(.if (•, M, DM)).ifp(', M, Dm) 

(4.3) + 

= ~ D^)]^ (d(-^p(l u, Dm)) - 

By mutual orthogonality of the projections in (2.1), the left and right hand side 
of (4.3) are normal to each other. Hence, they both vanish and we may split the 
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system to two components. We renormalise the second half of (4.3) by multiplying 
by TO — 1 to obtain 

D(.if (•, u, Du))^p{-, u, Du) 

(4.4) + ' [.Sfp(•,•», DM)]'^(D(.gp(-,M, Du)) -.if^(-,u, Du)^ =0, 

.if(-, u, Du)[.ifp(-, u, Du)]"*" ^D(jfp(-, M, Du)) — .if,,(-, u, Du)^ = 0. 

As TO 00 , we obtain the complete ODE system in L°° which after an expansion 
of derivatives and summation of the systems becomes (1.6)-(1.8). 

The degenerate elliptic case of the equations in L°°. Unfortunately, as we 
have already explained it is not in general possible to go much further without 
imposing the structural condition (1.13) on the hamiltonian .if. The problem is 
that the system fails to be degenerate elliptic in the sense needed for existence. 
In particular, the coefficient of (1.6)-(1.8) may be discontinuous at points where 
.ifp(-,u,Du) = 0. Assumption (1.13) forces the matrices an J^pp to com¬ 

mute and also makes the coefficients continuous by allowing them to match after a 
rescaling. 

We now formally derive (1.14) and also the Euler-Lagrange equations of the 
L™ functional (1-19) in the expanded form in which it will be used later, under 
the assumptions (1.15)-(1.16). There is no loss of generality in assuming > 1 
since if it is bounded below, we can always add a positive constant to and the 
equations remain the same because additive constants commute with the supremal 
functional (and this constant also regularises the minimisers of the respective L"* 
functional). In order to derive the equations, we first differentiate (1.13) and for the 
sake of brevity we suppress the argument {x, rj, P) of .if and (a;, 77, 4 |P — ^ [x, 77 )^) 
of and their respective derivatives: 


.ifp„ = 

■^Vc = - ^p{P-'^{x,r)))^i<yr,^{x,r]), 

= -^X - J^p{P - 'f'{x,r]))^iCfx{x,r]), 

■^Pa.Pf3= -^ppiP - ^{P - y(x,ri))^ + J^fpSap, 


■^P^x = -J^pPaxix,'!!) + {P - yix,!]))^ 


yPpx - J^piP -'P{x,v))P'rx{x,v) 


■^Pc-7/3 = -’^pPapfiix,!!) -k {P -y{x,l]))^ 



jepp{P - y{x,r]))^'rjpffix,r]) 


Recall now that (4.4) comprises the Euler-Lagrange equations of (1.9) in expanded 
form. Since by assumption Jifp > 0, we have the identities 


(4.5) 


[ji^p{P-y{x,r^))]^ = [P-r(cr,77)]^, 

[jy’p{P-y{x,r,))]^ = [P-r (37,77)]-". 


By grouping terms, setting 
(4.6) 


Wu := Du — y{-,u) 
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and omitting the argument (•, u, of and its derivatives, after a calcula¬ 

tion we have 




+ M’ppiWu) ® (Wu)^ + (J^pf(Wu) 0 (Wu) 


(4.7) .(^D\-rp(-,u)Du-%(-,u)^ + 




m — 1 


.^pCWu) ' 


-.Wp ^ {jffpp ■ Dm -h Jfpa:) Wuj + Jif’p{j^x + J^p -BujWu = 0. 
Similarly, by the identity (4.5) and since > 0, we have 

[Wu]^ (M’pl + J^pp{Wu)®(Wu)){p^u-'rp{-,u)V>u-%{-,u) 


(4.8) 


- Mp + ^p{Wuy-fpi-^u) + {.^pp-V)u + Mp^)Wu 


= 0. 


= 0 . 


Since the projection annihilates W (4.8) simplifies to 

(4.9) 

[ITm]-^ ^p(p^u-'fp{-,u)V>u-'f^{-,u'^ -M’p + M^p^WuY rp^-.u) 

We now observe that in view of the identities (2.1), the systems (4.7) and (4.9) can 
be matched and the mutually orthogonal coefficients \W'u\^ and [I^m]-*- add to the 
identity. By multiplying (4.9) by and adding it to (4.7), we obtain 


M’{M’p + ^pp\Wu\%Wu\ 


[^J^pf\Wu\^l 


V>{Wr 


(4.10) 


TO — 1 

-M’p^ M’p^Wuf 'fp{-,u) [M’pp-u^ M’p^){:Wu) 


M’[Wu 
TO — 1 

^ycx 


T 


^ -f ^p{^ 


-h J^p-V>u^Wu - .j>ifp\Wu\^[Wu]-^[.J>ifp - M’p{Wu)^yp{-,u)) = 0. 


The ODE system (4.10) is the Euler-Lagrange equation of the functional (1.19) in 
expanded form where for sake of brevity we have defined (4.6) and suppressed the 
dependence on the arguments (•, u, ^jWuj^) of and J^p, J^p, J^x- 

We also note that the coefficients which are of order 0(^^) are discontinuous, 
but this causes no problems since the terms involving these will be annihilated as 
TO —>■ oo. By letting to —>■ oo in (4.9) we obtain (1.6) with given by (1.14) and 
by (4.6). We finally rewrite the equations in a form which is more malleable 
for our proofs later. By setting 


(4.11) 


E“(-,m,Dm) :=- J^p(^J^x + J^p-Du^Wu 

+ {Jfpf\Wu\^[Wu]^(^J>fp - JfpiWu)^'rp{;u)'^, 
/“(•,m,Du) ■= - J^[Wu]^ (^ - J^p + M’p{Wu)^rp{-,u) 

+ {M’pp ■ u + ^px) Wv^, 


(4.12) 
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(4.13) A°°{-,u,Bu) := + J>fpp\Wu\^)[Wu]^, 

(4.10) can be written as 


(4.14) 

^“(•,m,Dm) 




u\^I 


B{Wu) = F^{-,u,Bu) 


TO — 1 
and (1.6) as 

(4.15) =F“(-,u,Du) 

where F°°,f°°,A°° are given by (4.11)-(4.13) and 'Wu by (4.6). 


TO — 1 


5. Existence of vectorial Absolute Minimisers 

In this section we establish item (1) of our main result Theorem 1.1 by proving 
existence of a mapping u°° G K^) which satisfies (1.2) for (1.17). 

Lemma 5.1 (Existence of minimisers and convergence). Let satisfy the 

assumptions of Theorem 1.1. Then, for any affine mapping h : K —> and any 
TO S N, the functional (1-19) has a minimiser u™ over the space ITjj^’^’”(r2,K^). 
Moreover, we have the estimate 

(5.1) ||w||vEi.2'"(n) < c(^Em{u,VL)^ + max|&| + 1^ 

for any u G M'^), where C > 0 depends only on the assumptions and the 

length of LI. In addition, there is a subsequence {mk)f° and u°° G W^’°° 
such that for any q> 1 

f u™- ^u°°, znC°(n,R^), 

\ Du™ -^ Du°°, in L«(D,K^), 

as TOfc —t oo, and also 

(5.2) ||R°°||ivi.=o(f2) < C. 

Finally, for any A C fl measurable with |A| > 0, we have the lower semicontinuity 
inequality 

(5.3) E^{u°°,A) < liminfE„(u™,A)^. 

m—>-oo 


Proof of Lemma 5.1. Step 1. We begin with some elementary inequalities we 
use in the sequel. For any t>0, 0<a<l and e > 0, Young’s inequality gives 

(5.4) t°‘ < et + ^ '*(1“®)- 

Moreover, for any P,V G R^ and 0 < 5 < 1, we also have 

(5.5) {1 - S)\P\^ < \P - V\^ + hv\^. 

d 

Finally, for any u G W^’^"‘(L2, R-^), we have the following Poincare inequality whose 
constant is uniform in m GN: 

||u||l 2™(0) < 2(|D| + l)(||Du||i2„(0) + max|u|y 


(5.6) 
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Indeed, in order to see (5.6), suppose u is smooth and since \u{x) — u{y)\ < |Du|, 

for y € dfl by Holder inequality we have 

^ ^ 2m 


/ iDul + maxlul ) 

Jq an ' 7 


< 2 


2m-l 


< (2(|f^| + l)) 


I |Dm| 
n 

2m-l 


+ max u 

on 


2m 


|Du| 


+ maxIuP"^ 

an 


which leads to (5.6). 

Step 2. We now show that the functional Em is weakly lower semicontinuous in 
Indeed, by setting 

/ li i2N™ 

(5.7) H(x,ri,P):=Jf’(^x,r],-\P-y(x,ri)\j , 

we have for the hessian with respect to P that (we suppress again the arguments 
of and its derivatives) 

Hpp{x,y,P) + (^M’J^pp + {m-l)[^pf^ 

.{P-r{;u))®{P-P{;u)) . 

By (1.15) and since the projection [P — y{-,u)Y satishes the matrix inequality 
[P — 7^(-,'«)]^< I, we obtain 

Hpp{x,r],P) > mjy"^-'^\jyjypl + jyjrpp{P-y{-,u))(S){P-y{-,u)) 

(5.8) > (M’pl + (co -M’p)[P- 

> m(^JfpI + {cq-^ p)[P 

> mcol. 


The conclusion now follows by standard lower semicontinuity results (e.g. [D, GM]). 
Step 3. Now we derive the energy estimate which guarantees the coercivity of Em- 
By our assumptions on there is a p S [0,p] such that 

jy{x,r],p) = Jifp{x,r],p)p + jy{x,p,0) > c^p + 1 . 

Hence, by using (5.5) the above gives 

(5.9) jf^{x,p,\\p-y{x,p)\') > |(i-5)|pp - 

Then, by (1.16) and (5.4)-(5.5), for cr > 0 small we have 

Jff’[x,p,l\p-y{x,p)\^) > |(i-<5)|Pp-^(1 + 77)2 

> |(1-<5)17^1" - - C(S,a,a), 

where C(S, a, a) denotes a constant depending only on the numbers <5, a, a. We now 
select 6 := 1/2, a := 2coe > 0 to find 

^(x,V,^\P -y{x,p)\^^ > j\P\^ - £\v\^ - C{e,a). 
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Hence, for any m S N by the Holder inequality and the above estimate, we have 

Consequently, for any u G K^), by integrating over H and by utilising 

(5.6) and (1.19), we deduce 

3(g)™ [ |Dw|2- < E^{u,n) +e^ ( |u|2- + C(e,a)2-|r!| 

t/ r2 r2 

< Em{u,9) + C(e,a)2-|H| 

+ c™(2(|H| + 1))^”^ {max|6p™ + ^ 

Hence, we have obtained the estimate 


{(S)™" (4(l^^l + lf^)™} / iD^^r < E^{u,n) + C2-(rnax|6r + l) 

i/ 

where C above depends on e,a, H. By choosing e := Co /(3 2 ®(|H| + 1 )^), we get 

and since limm_>oo (l — 2“"')^^™ = 1, the desired estimate (5.1) ensues. 

Step 4. We show existence of minimisers and convergence by using ideas of [BJWl]. 
We have the a priori energy bounds (recall the notation (4.6)) 


mi{Em{v,n)^ : < Em{b,n) 

^(•,6, |ir6|V2) 

|H|5k||^(.,5,|ir6|V2) 


L“(n) 


and Em{v,Vl) > 0, for any v G ’^'"(H, R-^). Hence, by standard results (see 
e.g. [D, GM]), there exists a minimiser u™ of the functional Em in R-^). 

Moreover, by (5.1) and (5.11) we have the bound 


(5.10) lh'"llvuu 2 m(Q) < C (^sup=^(-,6, |#^6p/2)'“ + inax|6| + 1^ . 

Let (7(0,6) denote the right hand side of (5.10). Then, for any q G [2,m], we have 


(5.11) 


||u"*||vUl,2,(n) < |0|29 2r,i 

< |0|5L-2kc'(f2^6). 


Hence, for any g > 1 fixed, the sequence ( 11 ™)“ is bounded in W^^’^'^(0, R^). As 

such, there exists u°° G 1 ’^«(0, R") satisfying u™ -^ u°° in W^^’^'^(0, R^) 

along a subsequence -G 00 . By letting m— 7>ooin(5.11) along the subsequence, 
the weak lower semicontinuity of the L^'^(0,R^) norm implies 

||w°°||vuu2<!(n) < |0| 251 (7(0,6). 

By letting now g —>■ 00 , we derive the desired bound for u°°. 
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Step 5. We finally show (5.3) by using ideas of [BJWl]. Fix A C n a measurable 

set with |A| > 0. By recalling that rt"* -^ as to —>■ oo along a subsequence in 

L'^{A,'R^) for any g > 1, by weak lower semicontinuity we have 

Eoo{u°°,A) = lim Eq{u°°,A)'i 

<J—>-oo 

< lim inf ("lim inf Fig(u™, A) 9 'j 

q—>oo V m—>oo J 

< liminf ("liminf 

q—¥oc) \ m—¥oo J 

< lim'mi Ern{u!^, A)^. 
m—xx) 

The lemma ensues. □ 


Now we are ready to establish the existence of Absolute Minimisers of the func¬ 
tional (1.17). 


Proof of (1) of Theorem 1.1. We continue from the proof of the previous lemma. 
Fix n' <E n. Since fl' is a countable disjoint union of intervals, it suffices to assume 
that n' = (a, b) d 11. We fix ^ S VFq(( a, b), M.^) and set := u°° + </>■ Hence, 
in order to show that u°° is an Absolute Minimiser of (1.17) over H it suffices to 
show that 

E^{u°°,{a,b)) < Eao{il^°°,{a,b)). 

Note also that u°°{a) = V'°°(a) and u°°{b) = 'ijj°°{b). We now fix 0 < 7 , (5 < {b—a)/S 
and define the following map: 






X G {a,a + j), 

X G (a + 7 , 6 - (5), 
X G {b — 5,6 ), 


where to S N U { 00 }. Then, we have g 

' -!/'“(» + 7 ) - (a) 

7 

ip°°{b-5)-u^{b) 

^ - 5 ' ^ ’ 


[{a,b),R^) and 

on (a, 0 + 7 ) 
on (a -I- 7 , 6 — S) 

on (b — S, b) 


Note now that 


(5.12) ^ in W^’°°((a, 6 ),K^), as to ^ 00 . 

Indeed, since obviously 'ip'^’'^’^ —)• 6 ),IR^), it suffices to note that 

for a.e. x G (o, b) 




\u°°{a) -u"^{a)\ |m°°( 6 )-m”"( 6 )| 

X(a,a+7) T \{b—5^b) 


< (- + i 1 \\u^-u^' 

7 0 


L“(n) 


7 


6 
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and hence 


L“(n) V7 5 ) 




(5.13) 


L“(n) \7 5^ 

as m —>■ 00 along a subsequence. Since for each m S N is a minimiser of (1.19) 
over Ibj’^"*(12, K^), by recalling that ^t {a, 6 }, minimality and Holder 

inequality give 

< £;„(V-™’^’^(a,6))- 

< (6-a)^^;oo(I/’™’^’^(a,6)). 

On the other hand, since 

Eoo (a, b)) = max jiloo (a, a + 7 )), 

i?oo(V'™’^’^(&-^, 6 ))} 

and = . 0 OO on (a + 7 , 6 — J), we have 

Eao{i^”"’'^'\ia,b)) < max|i;oo(i/’'"’'^’^,(a, a + 7 )), E^{tj}°°, {a,b)), 
(5.14) ^ 

i?oo(i/'™’^’^(5-<5,6))}. 

By combining (5.12)-(5.14) and (5.3) for A = {a,b), we get along a subsequence 
(rrii)^ that 

Eoo{u°°,{a,b)) < limmf (^max|£;oo(iA'"”'^’'^,(a,a + 7)), i;oo(i/’°°,(a, 5)), 

i?oo(V'™”'^’^(5-<5,6))}) 

which in turn gives 

Eoo{u°°,{a,b)) < iJia,x^E^{ip°°,{a,b)), E^{il;°°’'^’\{a,a + j)), 

(6-^,6))}). 

Further, since 

D^oo, 7.5 = on ( 0,0 + 7 ), 

D^oo. 7.5 ^ d1’-^i/>~( 6), on ( 6 -(5, 6 ), 


(5.15) 


we have 
(5.16) 


£^ 00 ( 0 , 0 + 7 )) = max 

[a.a+7] V / 

E^U°^’'^’^,{b-S,b)) = max 

[&—5,61 V 


By (5.15)-(5.16), we see that it suffices to show that there exist infinitesimal se¬ 
quences ( 7 i)^i and ((5i))^i such that 


(5.17) 


£'oo(V'°°) (o; 5)) >max|limsup max .if ("•,'0°°’'’'’’'^% D^’'''‘'0°°(a)'), 
1 i-s-oo la,a+ji] V / 

limsup max .if('•, 
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The rest of the proof is devoted to establishing (5.17) and this will complete the 
proof. 

We illustrate the idea of the proof of (5.17) by assuming first that exists 

on [a, b]. In this special case, we have 

[a,&] ' 

> max|.if^a,i/>'^(a),Di/;'^(a)^, .if ^5,'0°°(&),Di/!° 

Further, since the difference quotients satisfy 

^ Di/;~(a), Di/>~(6), 

as 7 , 5 —>■ 0, and also 

^00,7.5(2.) _^oo(^) -as 7 ^ 0 , 

-^0, as 5^0, 


(5.18) 


max 

a<ai<a+7 


max 

b—6<.x<.b 


we obtain 

lim max = Jffa,V'°°(a),DV^~(a)V 

7-10 [a,a+7] V /V / 

lim max ,B^’-^fj°°{b)] = .ifffe,i/>“(6),Di/>“(6)V 

5->-0 [&-i5,b] V /V / 

By putting these together we are led to (5.17). 

Now we return to the general case. Fix u G IF^’“(fI, K-^), x G [a, 6] and e > 0 
small and set 

Ag{x) := [x — e,x + s\(^[a,b\. 

Then, we claim that there is an increasing modulus of continuity lo G C^(0,oo) 
with w(0+) = 0 such that 

(5.19) E(u, As{x)) > ess sup .iffa:,M(a;),Du(?/)^ — uj{s). 

yeAe{x) ^ f 

In order to see (5.19), note that for a.e. y G A^(x) we have \x — y\ < e and hence by 
the continuity of .if and the essential boundedness of Du, there is an oj such that 


.if 


(^x,uix),I)uiy)^ - Af(^y,u{y),I)u{y)j < a;(e). 


for a.e. y G ^^(x). Hence, (5.19) ensues. Now we claim that 

(5.20) sup (limsup i —';F(.,'u) ) < ess sup ^ Dit — iF(., u) 

A,{x) V t^O 2 J A,{x) 2 

In order to see (5.20), it suffices to apply the inequality 

v{y + t) -v{y) 


t 


to the Lipschitz map 


<esssup|Dv|, y,y+ t G Ai:{x), t 0, 

Aeix) 


v{y) := u{y) - f y{t,u{t))dt 
J a 















GENERALISED SOLUTIONS AND VECTORIAL CALCULUS OF VARIATIONS IN L' 


19 


and note the identities 


nv{y) = Du(i/) - y{y,u{y)), a.e. y G Ae{x), 

1 /*y+i 

T)^'*v{y) = V>^'^u{y) - ^ J 'y{t,u{t))dt 

= D^’*'u(y) - 'r{y,u{y)) + o(l), as t 0, y G A^{x). 


Hence, (5.20) holds true. Now, we combine (1.13), (5.19) and (5.20) together with 
(1.15) to find that for any fixed u G K^), x G [a,b] and e > 0 small we 

have 

Eoo{u,{a,b)) > E^{u,As{x)) 

/ 1 2 
= ess sup y, u{y), - Du - r(y, u{y)) 

A,(x) \ 2 

> ess sup ix,u{x), ^|Du — {y,u{y))\^\ — uj(e) 

yeAe(x) \ 2 / 

/ 1 2\ 

= a;,u(a;),ess sup - Du(?/) — u(?/)) — a;(e) 


y^Ae(x) 

> M’ X, u{x), sup 

\ yeA^(x) 


limsup i|D^’*u(y) - 'E{y,u{y))\‘^ 


= sup 

y^Aeix) 


t —^0 


limsup ( x,u(x), -|D^’‘u(?/) — 'E(y,u(y))j 


- w(e) 


- w(e) 


> limsup ( ai,u(a;), i|D^’‘u(x) —';^(a;,u(x))|^ ) — uj{e) 

t-fO 


and by letting £ —>• 0 , we get 


(5.21) Eoo(u,{a,b)) > limsup ( x, u(a;), ^ |D^’‘u(a;) — l^(a;, u(a:))P] , 

t-j-o \ 2 ' ' / 


for any fixed u G 1H^’°°(H,IR'^) and x G [a, b]. Note now that since 

|d1’*u(x)| < ||Du||i=o(f 2 ), xG{a,b), t ^ 0, 

for any infinitesimal sequence (ti(ai))“i there is a subsequence denoted again by 
the same symbol such that 

(5.22) the limit lim D^’‘’*'“^u(a;) exists in K'^. 


By (5.21)-(5.22) and the continuity of we find that 


(5.23) 


Eaa{u,{a,b)) > limsup ( a;,u(a;),-|D^’*‘^'^^ u(x) — 1 ^( 2 ;,u(a;))| 

i—>oo \ ^ 

= Jif (x,u{x), ^ lim D^’*‘^'^^u(a;) — Y{x,u{x)) ) . 

y 2 i-^-oo J 


Now we apply (5.23) to 

u = X = a, b 

to infer that there exist sequences ( 71)^1 and (^i))^i such that 

(5.24) the limits lim {a) and lim (b) exist in 

i—^oo >-oo 
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and also 


Eooii ^'^> max ■( ( a,'i/’°°(a) i lim — l^(a,^°°(a)) 

' ' ’ ' 2 2—>-00 


(5.25) 


z l—¥00 


On the other hand, by (5.16), (5.18) and (5.24), for 7 = 7 i and S = 5i we have 

lim (a, a + 7 i)) = lim max .if (•, D^’'''‘'!/'°°(a)') 

WOO WOO [a,a+ 7 il V / 


> [a,a+7i] 

(5.26) =M’ (a,il;°^{a),\ lim Di’T'-^/>“(a) - r(a, V'°°(a))'^ 

Y 2 2—>^00 

and similarly 

lim = lim max 

2 —>-oo 2 —>-oo [6—(5^ ,6] V 

(5.27) =je (b,-il)°^{b),\ lim d1’-^-V^°°(6)-7(6,^/>°°(6))'^ 

Y 2 i—^00 

By putting together (5.25)-(5.27) we see that (5.17) ensues and so does item (1) of 
Theorem 1.1. □ 


6. Existence of P-solutions to the equations in L°° 

In this section we establish items (2)-(4) of Theorem 1.1. We begin by showing 
that the minimisers obtained in the previous section actually are weak solutions of 
the respective Euler-Lagrange equations. 

Lemma 6.1 (Weak solutions of the L™ equations). Let Jtf, 7, il, b satisfy the as¬ 
sumptions of Theorem 1.1 and let (u™)^ be the sequence of minimisers constructed 
in Lemma 5.1. Then, each u™ is a weak solution in Wf,’ ”^(12,K^) of the Euler- 
Lagrange equations of (1.19) on LI: 



where Wu is given by (4.6). 


Proof of Lemma 6.1. Let IT be given by (5.7). By suppressing for brevity the 
arguments of , J^p, we have 

Hp{x,T],P) = 

Hp{x,T],P) = - JtTp{P - 'r{x,r])Y'rp{x,rf)^ 

and the ODE system (6.1) can be written compactly as 
(6.2) D(i7p(-,u,Du)) = Hp{-,u,Du). 
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By (1.15), we have J^{x,ri,p) < (7(|7?|)(1 + p) and hence 

< C(h|)(l + |P-1 ^(t,77)|“). 

Further, by (1.15), (1.16) and the above, we easily obtain the bounds 

(6.3) \Hp{x,p,P)\ < C(|r?|)(l + 

(6.4) \H^ix,p,P)\ < CM){1 + |Pp™). 

By standard results (see e.g. [D]), (6.3)-(6.4) imply that the functional is Gateaux 
differentiable on M^) and the lemma follows. □ 

Now we show that the weak solutions u™ of the respective Euler-Lagrange equa¬ 
tions actually are smooth solutions. This will imply that the formal calculations of 
the previous section in the derivation of (4.10) make rigorous sense. 

Lemma 6.2 (C^ regularity). Let (it™)“ be the sequence of minimisers of Lemma 
6.2. Then, each u™ is a classical solution in (7^(17,K^) of the Euler-Lagrange 
equation (6.1), and hence of the expanded form (4.10) of the same equation. 


Proof of Lemma 6.2. Fix m > 2 and let us drop the superscripts and denote u"* 
by just u. The first step is to prove higher local integrability and then bound the 
difference quotients of Du in LjQj,(n, R-^). Let us fix g S N and C, S C^iLl) with 
0 < C < 1. By recalling (4.6), we set: 

(6.5) fiix) := C,{x) f f{t)\Wu{t)\‘^Wu{t)dt, x G fl. 

J inf fi 

Then, 4> G and 

D(j){x) = C'^{x)\Wu{x)\‘^Wu{x) -I-D()(a:) f f{t)\Wu{t)\‘^Wu{t)dt, 

J inf f2 

for a.e. a; G 17. Suppose now that q < 2m — 1. Then, since Du G L^™(17,R'^), 
we have that f) G W^^’^'"(17, K'^). By inserting the test function (p in the weak 
formulation of the system (6.2) (i.e. (6.1)) and by suppressing again the arguments 
for the sake of brevity, we have 





C^|iru|«iru -k DC / c,\WuWWt 
J inf n 


-k 



C [ C\^u\'^Wu 

. J inf r2 


0 . 


By (1.15)-(1.16), we have J^p > cq and 2h > co\Wu\'^. By using the bounds (6.3), 
(6.4) (where H is given by (5.7)), that C < 1 ^nd the elementary inequalities 


/. I/I ^ 

/ inf Q 
± 2 m—l 




in 

4.2m 


I/I, 


G 17, / G L^(17), 


< -k 1, 


t > 0 
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we have 



which gives 

/ <c(||w|u=o(n)) ( f ciWui^^+A. 

Jn \Jn J 

.J {|DC|(l + |#'up'"-i) + c(i + I^m|2-)|. 

Hence, we have obtained 

(6.6) f C"|irup-+« < C(||n|U=o(n)) ( [ [ l + liTwp-. 

Jn \Jn / Jn 

In view (6.6), by taking g + 1 = 2m we have 'Wu G L]^™“^(r2,K^). Hence, we can 
iterate and choose g + 1 = 4m — 1 to find that (j) £ M'^) which makes it 

admissible and we can repeat the process. Hence, by applying the estimate again 
we infer that WuG By induction, the estimate holds for all integers 

of the form 

q = (2m — l)fc, fc £ N 

and we obtain that Wu G n^i view of (4.6) and since u G 

C°(H,IR'^), we conclude that 


r=l 


The next step is to prove that D^’‘Du is bounded in The idea is classical, but 
we provide the arguments for the sake of completeness. To this end, we test in the 
weak formulation of (6.2) against difference quotients of the form 

J, T^i,—i 'v ^ ^ ^ rt(x -f t) li(x) , / r, 

0 := -D ’ (C D ’ uj , c G (H), D u{x) = - - -, t ^ 0. 

Let H be given by (5.7). Then, for £ > 0 and t small, we have 


I ■■= J Di>*(Hp(-,u,Du)) • (c^Di’‘Du + 2CDCD1 Gu) 

(6.7) < ^|i7„(-,u,Du)|(c|DiGu| + 

<k( f C\Hrj{;U,Bu)\^ + f ClDup) +e [ C^IDI’^DuI^. 

\Jq Jn J Jn 


for some constant K > 0 independent of t. By using the inequality Hpp > cqI and 
the identity 
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+ Xu{x + i) + (1 — A)u(a;), ADu(a; + t) + (1 — A)Du(a;)^D^’‘M(a;) 

+ Xu{x + <) + (! — X)u{x), XDu{x + t) + (1 — A)Du(a;)^ j- dX 

(where Hp^^ ^ denotes difference quotient with respect to the x variable), we have 

(6.8) i - C'(||u|U»(a)) f C|P(|Du|)| 

^ Jq Jn 

where K > 0 is a constant independent of t, whilst P is a polynomial expression and 
it is a consequence of (1.15)-(1.16). Since u G and Dm G K-^) for 

all r > 1, by (6.7) and (6.8) we obtain that u G bPjQ’^(D,K^). Thus, the calculations 
in the derivation of the expanded form of the system make sense a.e. on Q. Since 
Hpp is a strictly positive matrix, by a standard bootstrap argument in the system 
we obtain that u G and the lemma follows. □ 


Now we may prove the remaining assertions of our main result. 


Proof of items (2)-(4) of Theorem 1.1. In view of Lemmas 5.1, 6.1, 6.2, 
let (m"*)J° denote the sequence of minimisers in C'°(D,IR'^) n (7^(12,IR-^) of the 
functionals (1.19) over the spaces M^). Then, along a subsequence 


(6.9) 


M 

Dm”’ 


^ Dm°°, 


in C0(D,K^), 
in L'J(D,]R^), for all <? > 1, 


as m —>■ 00 , and the limit satisfies m°° G ’°°(D,]R^). Moreover, each m™ is a 
classical solution of the system (4.10), or equivalently of (4.14) with f°°,F°°,A°° 
given by (4.11)-(4.13). The goal is to show that the limit map u°° is a D-solution 
of the system (1.6) with Fao given by (1.14) (or equivalently (4.15)) and also that 
M°° = 6 on dfl. We begin by observing that the boundary condition is satisfied 
as a result of the uniform convergence on D. Moreover, by recalling (4.6) and by 
multiplying (4.14) with D(1 ^m"’), we obtain 


(•,M, ^ I /Tm"’II ir m”’I"/| : B{Wu^) 


( 6 . 10 ) 


m — 
®D(>rM"’) = 

< 


/°°(-,m"’,Dm”’) 
m — 1 

/°°(-,m"’,Dm”’) 


+ F“(-,m“,Dm"’) ) •D(irM”’) 


m — 1 


+ F~(-,m”’,Dm’") 


\D{Wi 


By (1.15) we have JFp > cq. In addition, by (4.13) the matrix map A°° i 
negative. Hence (6.10) gives the estimate 

(6.11) IiTm’”! D(irM”’) < ^ ■' 1 +F“(-,m”’,Dm"’) 

By using the elementary inequality 

|D(|/|3)| < 3||/pD/|, /gC1(D,M^), 

(6.11) gives the estimate 


is non- 


( 6 . 12 ) 


D 


(IiTm 


jn |3\ 


- ^2 


fOQ(, m r).ym\ 

(•, m"* , Dm™ ) 

m — 1 
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By (6.12), (6.9) and the form of the right hand side given by (4.11), (4.12), we have 
that the sequence 

^0 X3) ^ m. IR I-i->. m. m. \ 1 3 


:= \Wu^Y = pw™-r(-,n'") 


is bounded in for any q>l. Hence, by the compactness of the imbedding 

(E C^(yi)^ there is a continuous non-negative function such that 


in C7°(fl), 


V —> V 

along perhaps a further subsequence as m —)■ oo. We claim that 


(6.14) 




,oo |3 


= - r(-,u°°)| < v°°, a.e. on n. 


Indeed, by (6.9) and the weak lower semi-continuity of the norm, for every x G 
and r > 0 fixed we have that 


I 




\Wu 


(6.15) 


< lim inf — 

m—^oo 2r 


lim — 

m—^oc) ZT 


/*tc+r 


\Wu^ 


x—r 

nx-\-r 




/'X-\-r 


By passing to the limit as r —)■ 0 in (6.15), the Lebesgue differentiation theorem 
implies that the inequality (6.14) is valid a.e. on 12. We now set 

:= {x Gil : v°°{x) > O}. 

By the continuity of v°°, il°° is open in H, the set 12 \ 12°° is closed in 12 and 
f2\f2°° = {x G 12 : v°°{x) = 0}. 

By (6.14), we have 


(6.16) 


I Wu°° 1=0, a.e. on 12 \ 12°°. 


On the other hand, since n™ —> v°° in (7° (12), for any U (s f2°°, there is a cto > 0 
and an m{U) G N such that for all m > m{U), we have w™ > ctq on t/ and hence 
by (6.13) 

(6.17) > (o-o)*, onU. 

By (6.17), (6.15) and (6.12), we have 


(6.18) D(iru” 


< 


(co)^(o-o)3 


/°°(-,m”*,Du™) 


m — 1 

\2^,m 


-FF°°(-,u™,Dw'") 


on 12'. 


By (6.18) and (6.9) we have that is bounded in ,'R^). Hence, we 

have that 



- M°°, 

in C'°(D°°,K^), 


Du’" 

—^ Du°°, 

in L«„^(D°°,K^), 

for all q > 1, 

D2ym 

-- D2u°°, 

in 

for all (j > 1. 


Thus, by passing to the limit in the ODE system (4.10) as m —>■ oo along a sub¬ 
sequence, we have that the restriction of u°° over the open set 0°° is a strong a.e. 
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solution of (4.15) on 11°°. By bootstrapping in the equation, we have that actually 
u°° G (fl°°, R-^). On the other hand, we have that 

\Wu°°\=0, a.e. onO\0 °°. 

Hence, if the set H \ has non-trivial topological interior, by differentiating the 
relation we have that exists a.e. on the interior of the open 

set O \ and by bootstrapping again we see that u°° G (7^(int(0 \ 

Putting the above together, we have that exists and is continuous on the 

open set Ooo defined in the statement of the theorem which is the union of 0°° and 
of the interior of fl \ 0“: 

G Ooo = 0~uint(0\0°°). 

We now show that is a P-solution of (4.15) on H (Definitions 3.1-3.2). Let 
be the first difference quotients of Du°° along a sequence —>■ 0 as 

i —>■ oo and let 'D^u°° be a diffuse 2nd derivative of u°° arising from the subsequential 
weak* convergence of the difference quotients, that is 

in<r(0 ,l^), 

as j —^ oo, in the space of Young measures from H C K into the 1-point compact- 
ification R^ = R^ U {oo}. By the regularity of u°° on and Lemma ??, the 
restriction of any diffuse 2nd derivative on is the Dirac mass at the second 
derivatives: 


(6.19) 'D'^u°°{x) = (5d2„cx=(2.), for a.e. x G 


Hence, u°° is I?-solution on 12°°, since it is a strong solution on this subdomain. 
Consequently, for a.e. x G 12°° C 12 and any X G supp* (l22u°°(a;)) we have 


( 6 . 20 ) 


(x, (x), ^ I Wu'^ (x) I I (x) 

= i^°°(x,u°°(x),Du“(x)). 


X-D(r(-,u“))(x) 


Thus, u°° is a I?-solution of (1.6) with Too given by (1.14) (i.e. (4.15)). On the 
other hand, since |1 ^m°°| = 0, a.e. on 12 \ 12°°, for a.e. x G 12 \ 12°° and any 
X G supp.„(l2^u°°(x)) we have we have 


( 6 . 21 ) 


<1 


X, u 


\x),-\Wu° 


(x)n|iru°°(x)|^ Y-D(ir(-,M°°)(x) 


= 0 . 


Also, by (4.11) we see that the right hand side of (4.15) essentially vanishes on 
12 \ 12°° as well: 

(6.22) i^°°(-,u°°,Du°°) = 0, a.e. on 12 \ 12°°. 


By putting (6.20), (6.21), (6.22) together, we conclude that u°° is indeed a 12- 
solution of the Dirichlet problem for the fundamental equations in L°°, which is also 
a weak sequential limit of minimisers of the respective functionals as m —> oo 
in the W^'‘^ topology for any g > 1. In order to conclude it remains to establish 
the strong convergence of the derivatives of the sequence of minimisers to u°°. 
On the open set 12°° we have Du^ — Du°° in L^^^(12°°,R-^) and hence up to a 
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further subsequence we have Du'^{x) —>■ Du°°{x) for a.e. x G as m —>■ oo. On 
the closed set O \ , we have 



for any g G N. Since 

as m —>■ oo in (Q \ and also u™ —>■ u°° in C°(r 2 ,K^), we have that 

Du™ —S' Du°° in L^(fl \ ,M.^) along a subsequence as m —>■ oo. Conclusively, 

DvJ^{x) —S- Du°°{x), for a.e. a: G D as to —>■ oo along a sequence 
and also 

||Du™|U,(n) < C(g), 9 G N. 

Hence, if i? C O is measurable, we have the equi-integrability estimate 
I|Du™|U,(£) < ||Du™|U,+i(E)|i?|™ < C{q+l)\E\^. 

The conclusion of strong convergence now follows from the above and the Vital! 
convergence theorem. Theorem 1.1 has been established. □ 

Acknowledgement. The author is indebted to Jochen Brocker for the scientific 
discussions on the subject of variational Data Assimilation. 
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